Superconductivity in lithium, potassium and aluminium under extreme pressure: 

A first-principles study 
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Extreme pressure strongly affects the superconducting properties of "simple" elemental metals, 
like Li, K and Al. Pressure induces superconductivity in Li (as high as 17 K), while suppressing 
it in Al. We report first-principles investigations of the superconducting properties of dense Li, K 
and Al based on a recently proposed, parameter-free, method. Our results show an unprecedented 
agreement with experiments, assess the predictive power of the method over a wide range of densities 
and electron-phonon couplings, and provide predictions for K, where no experiments exist so far. 
More importantly, our results help uncovering the physics of the different behaviors of Li and Al in 
terms of phonon softening and Fermi surface nesting in Li. 

PACS numbers: 74.25.Jb, 74.25.Kc, 74.20.-z, 74.70. Ad, 71. 15. Mb, 74.62.Fj 



The effect of high pressure on phonon mediated su- 
perconductors has been the subject of many investiga- 
tions. These studies revealed a strong material depen- 
dence: while applied pressure suppresses superconduc- 
tivity in some materials, it favors it in others 0. Even 
in simple metals, the physics underlying pressure effects 
on the superconducting properties can be very compli- 
cated. For example, Li Q, 0, H H H H Q, [lfl 03 and 

ai m e h Il5| . behave in many circumstances like 

nearly free-electron gases, but they exhibit very different 
behaviors under pressure, still only partially understood 
within the Eliashberg theory [l6| . At ambient pressure, 
Al is a superconductor with T c = 1.18 K |l2j. Pres- 
sure rapidly reduces T c bringing it down to 0.075 K at 
6.2 GPa [L2j. Lithium, on the other hand, is a rather 
complex material: below 77 K and at zero pressure, it 
shows a martensitic transition to energetically compet- 
ing closed packed structures from 7.5 to 70 GPa it 
undergoes several structural transitions which suggest 
the presence of strong electron-phonon (e-ph) interac- 
tions. No sign of a transition to a superconducting state 
above 4 K was found up to rj 20 GPa while, at higher 
pressures, Li becomes a superconductor J3|, |4j, |5j, |6J . In 
the range 20-38.3 GPa, where Li crystallizes in an fee 
structure, experiments by Shimizu Q, Struzhkin [jj, and 
Deemyad |6J found that T c increases rapidly with pres- 
sure, reaching values around 12-17 K (the highest T c 
observed so far in any elemental superconductor). Also 
K undergoes several phase transitions, and is stable in 
the fee phase between 11.6 and 23 GPa|l8|. 

The pressure dependence of T c for Al, Li and K has 
been calculated by Dacorogna et al. |L4^. by Christensen 
and Novikov 0, HJ , and by Shi et al. respectively. 
For Al, Dacorogna et al. [l4| obtained a nearly satisfac- 
tory agreement with experiments [l2l Il3| . In the case of 
Li, Christensen et al. |?J used a rigid-muffin-tin approx- 



imation for the e-ph coupling constant A and /i*—0.13. 
Due to the empirical scaling of phonon frequencies y|, 
they obtained a much too high T c (45-75 K) unless the 
"non-standard" value of /z*ss 0.22 or an additional term 
modeling spin-fluctuations was used. For K, Shi et al. || 
obtained T c = 9 K at 13.5 GPa, using /i*=0.13. 

Note that in these reports, as in most other Eliashberg- 
based calculations, the electron-electron (e-e) repul- 
sion was treated semi-empirically through the Morel- 
Anderson pseudopotential /i*. The validity of this pro- 
cedure for the case of low density solids like Li has been 
questioned by Richardson and Ashcroft . Building on 
the seminal work of Oliveira, Gross and Kohn [2(J, and 
on further developments |2l| , some of us recently intro- 
duced an alternative approach to Eliashberg theory: an 
extension of the density functional theory to the super- 
conducting state (SCDFT) 22]. This theory is fully ab 
initio, and is capable of describing correctly the super- 
conducting properties of several elements |23j and com- 
pounds |24j. In the present paper we explore this very 
promising method to study the superconducting proper- 
ties of Li, K and Al under pressure. Furthermore, we 
provide a detailed description of the subtlety of super- 
conductivity in Li, where the incipient phase transitions 
produce a phonon softening and a very strong electron- 
phonon coupling, thus enhancing T c up to values unusu- 
ally large for simple elemental metals. Our results for 
Li and K confirm that a full treatment of electronic and 
phononic energy scales is required, in agreement with pre- 
vious arguments 0. We predict K to be superconduct- 
ing, with a T c up to sa 2 K in the experimental stability 
range of the fee structure, and up to w 11 K in the range 
of stable calculated phonon frequencies. 

Ground-state calculations were performed using the 
pseudopotential based code PWSCF within the lo- 
cal density approximation (LDA) to the density func- 
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FIG. 2: Upper panel: phonon dispersion of Li along the X — 
K — r line, at different pressures, for the lower frequency mode 
(frequencies below the zero axis denote imaginary values). 
Lower panel: electron- phonon coupling A 9i i and phonon line- 
width 7q,l. 



FIG. 1: Comparison between calculated and experimental 
critical temperatures for fcc-Al (upper panel), K (middle 
panel) and fcc-Li (lower panel). Numbers in parenthesis after 
'McMillan' indicate the /i* value. Different symbols with the 
same color refer to the same experimental report using dif- 
ferent setups. Vertical dashed lines indicate the experimental 
structural transition pressures for Li and K. The insets depict 
A vs pressure in GPa. 

tional theory. The validity of the pseudopotential ap- 
proach at high pressures was verified by comparison with 
all-electron methods. Phonon frequencies and e-ph cou- 
plings were obtained from density functional perturba- 
tion theory. The electron-phonon coupling spectral func- 
tion a 2 F{uj) and q-dependent phonon lincwidth were 
evaluated through a careful integration over the Fermi 
surface. The implementation of SCDFT has been re- 
ported elsewhere |22j, |23| . All systems were considered in 
the fee structure in their experimental ranges of stability. 

The LDA underestimates the equilibrium volume of Li 
and Al, as it can be observed by calculating the equation 
of state P(V) using Murnaghan's formula. To compen- 
sate for this systematic error, we apply a positive pressure 
shift of about 3.5 and 2 GPa for Al and Li, respectively. 
This is the amount required to match the experimental 
equation of state @, 12a, |23 • Note that these shifts will 
always be included in the values given in the following. 
No shift is necessary for K. 

In Fig. ^ we compare the calculated pressure depen- 
dence of T c for Al, K and Li with available experimental 
results. For Al, SCDFT calculations match exactly the 
experimental zero pressure T c =1.18 K, and reproduce 
the rapid decrease of the transition temperature. In the 
same figure we report the estimation of T c by means of 



McMillan's formula (using /i*=0.13, in agreement with 
previous studies [l4|L 

In the case of Li (lower panel of Fig. [TJ, despite the 
poor agreement among the four sets of experimental data, 
the most recent experiments 0, agree in: (i) Li is not 
superconducting at ambient pressure; (ii) T c is lower than 
4 K up to 20 GPa; (hi) T c then increases with pressure, 
reaching 14 K at about 30 GPa. The only exception to 
this behavior is the early report by Lin and Dunn 
In this pressure range, the results of Struzhkin [j| and 
of Deemyad Q are quite similar, while Shimizu et al. 
find lower values of T c . At even higher pressures, ex- 
periments show a quite complex behavior (see below). 
Within this scenario, our calculated SCDFT results are 
in excellent quantitative agreement with the most recent 
experiments yjlfj up to about 30 GPa. We find that Li is 
not superconducting up to 8 GPa, and that T c shows two 
different trends with pressure, a first region (8-20 GPa) 
in which T c increases at a rate of ~ 0.3 K/GPa, and a 
second region (20-30 GPa) at - 1.3 K/GPa. 

K shows a behavior quite similar to Li: beyond a pres- 
sure threshold (20 GPa) T c rises rapidly. In the range 
where phonons were found to be stable, it reaches « 11 
K at 29 GPa; the experimentally observed instability of 
the fee phase, however, limits this value to « 2 K at 
23 GPa. 

The differences between Al and Li can be understood 
by looking at the e-ph coupling as a function of pressure. 
In Al the phonon frequencies increase as the pressure rises 
(this corresponds to the normal stiffening of phonons 
with increased pressure). In addition, the Eliashbcrg 
spectral function a 2 F(uj) shows a high-frequency peak 
whose height decreases as a function of pressure. These 
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FIG. 3: Panel (a): 3D view of the Fermi surface of Li, at 
28.6 GPa, with a color scale indicating the value of A*,. The 
red color indicates high values of A. Panel (b): the FS cut 
on the plane of Panel (a), including the £-point and paral- 
lel to the (110) plane brings into evidence the FS nesting. 
The arrows represent nesting vectors. Panel (c): isosurface of 



A,;, 



= i = 5. F is at the vertices and at the center of the cube, 



X at the centers of the faces. 



factors contribute to the decrease of A (see the inset of 
Fig^) and consequently of the critical temperature T c . 

In Li, on the other hand, the phonon frequencies ex- 
hibit a quite different behavior. In Fig. [21 we present 
the phonon dispersion of the lowest branch along the 
X—K—T line of the Brillouin zone (BZ). For pressures up 
to 8 GPa, there is an increase of the phonon frequencies, 
as in the case of Al. However, as the pressure is raised fur- 
ther, the phonons near the if -point start to soften. The 
softening continues up to 33 GPa when this frequency 
becomes imaginary. A closer inspection of Fig. [5] reveals 
that already at around 30 GPa a phonon mode close to 
the T-point develops an imaginary frequency. We be- 
lieve that this marks the transition to the hRl phase, 
but further analysis would be required to fully validate 
this assumption. Although diffraction experiments per- 
formed at 180 K set the structural phase transition 
at 39 GPa (up to 42 GPa), it was observed recently that 
T c has a maximum at 30-33 GPa and drops drastically 
beyond that pressure 0, [||. This latter finding, so far 
unexplained, is consistent with our theoretical prediction 



of a complete phonon softening at around 30 GPa. 

In Fig. |2 we plot for the lowest frequency phonon 
the linewidth 7 g „ [the Fermi surface (FS) average of 
the e-ph coupling matrix elements] and the correspond- 
ing q— dependent electron-phonon coupling \ q ,„, at P = 
28.6 GPa. The quantity jq >u shows a peak close to K, and 
a broad maximum between K and X. This peak suggests 
the presence of FS nestings. To demonstrate this idea, we 
plot in Fig. |3| a cut of the FS parallel to the (110) plane 
which includes the L point. We clearly recognize nesting 
vectors (indicated by arrows) connecting fairly flat and 
parallel lines. The effects of nesting are remarkably en- 
hanced by the strong e-ph matrix elements. In Fig.^c), 
we show the isosurface A g ,„ = i = 5. As it can be seen, 
the extremely high-coupling regions form tubular struc- 
tures oriented parallel to the Cartesian axes, and cen- 
tered around q-values along the (110) direction, match- 
ing the FS nesting vectors indicated above. The observed 
phonon softening as a function of pressure is a direct 
consequence of the presence of the FS nesting. In turn, 
the progressive FS nesting with the increase of pressure 
is a consequence of the FS topological transition from 
a spherical free electron-like to a distorted anisotropic 
shape (see Fig. 01 . This topological transition is a mani- 
festation of an "s — p" transition of the electronic states 
near E F 0,0, E3] In particular, we can see in Fig.^a) 
that the k-resolved A*, is maximum on the rings of the 
FS around the L, arising from mostly p— like and strongly 
covalently bonded states. An analogous "s — <F transi- 
tion occurs in K. Also the electron localization function 
(ELF), which progressively increases with pressure, indi- 
cates the departure from the free-electron picture towards 
more "covalent" phases. This produces strong electron- 
phonon coupling that leads to the symmetry-breaking 
phase transitions at high pressure 0] and to the "paired 
bonded" structures close to 100 GPa pjof. 

Our calculation of a 2 F(uj) for Li leads to a value of 
A = 0.38. Using this number and the standard value 
^*=0.13 inside McMillan's formula we obtain T c =0.25 K 
at zero pressure, which is lower than the previously re- 
ported one , but in agreement with more recent calcu- 
lations 0- However, it is in complete disagreement with 
our SCDFT results and with experiments, both giving 
T c —0 K. In order to obtain this latter result from McMil- 
lan's equation we need to use fi* « 0.22. This value also 
describes quite well the behavior of T c with pressure. The 
same physical result was found for K (Fig. [Q, where a 
very similar value fi* ~ 0.23, has to be used to bring 
McMillan and SCDFT results into agreement for all pres- 
sures 33]. According to Richardson and Ashcroft [l9| . 
a large value of \x* is justified by the fact that the e- 
e interaction becomes unusually large at such densities. 
These authors pointed out that, in this density range, a 
full treatment of electrons and ions on the same foo ting 
is required. This is exactly what our method [22I Ejf 
achieves: The different energy scales of Coulomb repul- 
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sion and phonon-mediated attraction are fully included 
without any ad-hoc modeling. 

In order to further investigate the effect of the Coulomb 
repulsion, we solved the SCDFT gap equation in Li using 
for the e-e repulsion both the actual matrix elements of 
the Thomas-Fermi screened potential and the approxi- 
mate functional /C TF_SK introduced in Ref. |2^|. While 
this latter approximation failed for the localized a or- 
bitals in MgB2, it works very well for Li over the pressure 
range considered. Thus, while the strong e-ph coupling 
and the ELF indicate that covalency is becoming pro- 
gressively more important, the Coulomb repulsion can 
still be modeled as if electrons were delocalized. Then, 
why n* has to be so large? We can try to answer this 
question using SCDFT. By making a comparative anal- 
ysis between Al and Li we found that the FS average of 
the Coulomb potential multiplied by N(Ep) (which can 
be directly related to fi) is « 40% larger in Li than in 
Al, mostly because of the Van Hove peak in the den- 
sity of states of Li. This effect is present even after 
taking into account the renormalization of the Coulomb 
repulsion due to retardation effects. Thus, our results 
clearly demonstrate a larger value of the Coulomb repul- 
sion in Li than in Al. Since our functional, based on the 
Thomas-Fermi screened Coulomb repulsion, does not in- 
clude spin-fluctuations, we can clearly rule out that spin- 
fluctuations are required to explain the large Coulomb 
contributions and the anomalously high value of jj* . Fur- 
ther investigations at T = indicate that the gap func- 
tion is very sensitive to changes of the density of states at 
all energies, once again showing the difficulty to reduce 
the complexity of the Coulomb interaction into a single 
numerical parameter. 

In summary, the recently introduced SCDFT method 
allowed us to calculate the superconducting transition 
temperature of Al, K and Li under high pressure from 
first principles. The results obtained for Al and Li are in 
very good agreement with experiment, and account for 
the opposite behavior of these two metals under pressure. 
Furthermore, the increase of T c with pressure in Li is 
explained in terms of the strong e-ph coupling, which 
is due to changes in the topology of the Fermi surface, 
and is responsible for the observed structural instability. 
Finally, our results for K provide predictions intriguing 
enough to suggest experimental work on this system. 

apei 

)honon softening |jnj and 
in Li under pressure. 



Note: during the reviewing process of this paper two 
reports appeared confirming 
Fermi surface deformation 132 
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